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ABSTRACT 

Context. The Rosetta mission of the European Space Agency has been orbiting the comet 67P/Churyumov-Gerasimenko (67P) since 
August 2014 and is now in its escort phase. A large complement of scientific experiments designed to complete the most detailed 
study of a comet ever attempted are onboard Rosetta. 

Aims. We present results for the photometric and spectrophotometric properties of the nucleus of 67P derived from the OSIRIS 
imaging system, which consists of a Wide Angle Camera (WAC) and a Narrow Angle Camera (NAC). The observations presented 
here were performed during July and the beginning of August 2014, during the approach phase, when OSIRIS was mapping the 
surface of the comet with several filters at different phase angles (1.3°-54°). The resolution reached up to 2.1 m/px. 

Methods. The OSIRIS images were processed with the OSIRIS standard pipeline, then converted into I/F radiance factors and 
corrected for the illumination conditions at each pixel using the Lommel-Seeliger disk law. Color cubes of the surface were produced 
by stacking registered and illumination-corrected images. Furthermore, photometric analysis was performed both on disk-averaged 
photometry in several filters and on disk-resolved images acquired with the NAC orange filter, centered at 649 nm, using Hapke 
modeling. 

Results. The disk-averaged phase function of the nucleus of 67P shows a strong opposition surge with a G parameter value of - 
0.13+0.01 in the HG system formalism and an absolute magnitude bO) = 15.74+0.02 mag. The integrated spectrophotometry 
in 20 filters covering the 250-1000 nm wavelength range shows a red spectral behavior, without clear absorption bands except for 
a potential absorption centered at ~ 290 nm that is possibly due to SO 2 ice. The nucleus shows strong phase reddening, with disk- 
averaged spectral slopes increasing from 11%/(100 nm) to 16%/(100 nm) in the 1.3°-54° phase angle range. The geometric albedo 
of the comet is 6.5+0.2% at 649 nm, with local variations of up to ~ 16% in the Hapi region. From the disk-resolved images we 
computed the spectral slope together with local spectrophotometry and identified three distinct groups of regions (blue, moderately 
red, and red). The Hapi region is the brightest, the bluest in term of spectral slope, and the most active surface on the comet. Local 
spectrophotometry shows an enhancement of the flux in the 700-750 nm that is associated with coma emissions. 

Key words. Comets: individual: 67P/Churyumov-Gerasimenko, Methods: data analysis. Techniques: photometric 


1. Introduction 

The Rosetta spacecraft arrived on 6 August 2014 at comet 
67P/Churyumov-Gerasimenko (67P) after ten years of interplan¬ 
etary journey. The nucleus has been mapped by the Optical, 
Spectroscopic, and Infrared Remote Imaging System (OSIRIS) 
(Keller et al., 2007), which comprises a Narrow Angle Camera 
(NAC) for nucleus surface and dust studies, and a Wide Angle 
Camera (WAC) for the wide-field coma investigations. From the 
beginning of July, images obtained with different filters (from 
240 to 1000 nm) with the NAC camera allowed us to obtain 
color mapping of the comet nucleus with an increasing spatial 
resolution during the global mapping phase at the comet. Obser¬ 
vations will continue during the nominal mission until the end of 
December 2015. 


The very first results on the comet images and spectroscopy ob¬ 
tained with the OSIRIS and VIRTIS instruments (Sierks et al., 
2015; Thomas et al., 2015; Capaccioni et al., 2015) reveal that 
the nucleus of 67P has surface characteristics that are very differ¬ 
ent, in terms of shape, complex morphology, and spectral proper¬ 
ties, from those of the other cometary nuclei that were visited by 
space missions. Before the ESA Rosetta interplanetary mission, 
exploration of comets was considered by different space agen¬ 
cies. These began in 1986 with the flyby of comet IP/ Halley 
(Keller et al., 1986) by the ESA Ciotto mission, launched in 1985 
and followed by a flotilla of spacecraft from the Japanese and 
Russian space agencies. The Giotto Extended Mission (GEM) 
made the second flyby of comet 26P/Grigg-Skjellerup in July 
1992. In 2001, the NASA Deep Space 1 mission flew by comet 


Article number, page 1 ofm 



A&A proofs: manuscript no. 25901_xarv 


19P/Borrelly (Soderblom et al., 2002). The NASA Discovery 
program selected the missions Stardust and Deep Impact. The 
sample-return mission Stardust had flown by comet 81 P/Wild 2 
in 2004 (Brownlee et al., 2004), while in July 2005, the mission 
Deep Impact could observe the collision induced by an impactor 
on the nucleus of comet 9P/Tempel 1 (A Hearn et al, 2005). 
The comet 103P/Hartley 2 was visited in November 2010 by 
the redirected Deep Impact mission, named EPOXI (A Hearn et 
al., 2011), while Stardust, renamed NExT, was redirected to fly 
by comet 9P/Tempel 1 in Eebruary 2011 (Veverka et al., 2013). 
These missions (see Barucci et al. (2011) for a complete review) 
yielded an incredible knowledge of comets, unveiling their sur¬ 
face structures and many other properties. The known cometary 
nuclei are irregularly shaped, and each of them presents a vari¬ 
ety of different morphologies (e.g., depressions, ridges, chains 
of hills, smooth areas, rough terrains, layers, and craters) and 
different spectral properties. The Rosetta mission is the first mis¬ 
sion following and orbiting a comet from ~ 4 AU inbound to 2 
AU outbound including the perihelion passage at 1.24 AU. The 
OSIRIS imaging system is the first instrument with the capabil¬ 
ity of mapping a comet surface at such a high resolution (lower 
than 20 cm/px at best) with 20 Alters and covering various phase 
angle and illumination conditions. 

We here present results for the spectrophotometric properties 
of the nucleus of 67P derived from the OSIRIS observations ob¬ 
tained from July to mid-August 2014, during the comet approach 
phase and the first bound orbits. We present the global and local 
analysis of the photometric properties of the nucleus and spec¬ 
trophotometry. These data give indications for the properties of 
the comet nucleus and allow us to investigate the heterogene¬ 
ity of the nucleus at several scales, both in terms of albedo and 
composition. 

2. Observations and data reduction 

The NAC has a held of view of 2.2° x 2.2° and was designed to 
obtain high-resolution images with different Alters in the near- 
UV to near-IR range (240-1000 nm). These Alters are optimized 
for mineralogical studies of the nucleus. The WAC has a wide 
held of view (11.6° x 12.1°) and a set of narrow-band Alters de¬ 
voted to study gaseous species in the coma. The NAC and WAC 
cameras have unobstructed mirror systems with focal lengths of 
72 cm and 13 cm, respectively. Both cameras are equipped with 
2048x2048 pixel CCD detectors with a pixel size of 13.5 jum. 
The image scale is 3.9 "/px for the NAC and 20.5 "/px for the 
WAC. We refer to Keller et al. (2007) for a more detailed de¬ 
scription of the OSIRIS cameras. 

Erom mid-July onward, the comet was resolved, and OSIRIS 
has mapped the nucleus surface with increasing spatial resolu¬ 
tion. We present the results of OSIRIS observations obtained 
from the very first resolved images up to the beginning of Au¬ 
gust, when the spacecraft began to orbit the comet at a distance 
of 110 km. In this time frame, OSIRIS has mapped the sur¬ 
face of the comet with several Alters at different phase angles 
(1.3°-54°) and with a resolution of 2.1 m/px (Table[T]). The data 
were reduced using the OSIRIS standard pipeline up to level 

3, following the data reduction steps described in Kuppers et 
al. (2007) and Tubiana et al. (2015a), which include correction 
for bias, flat field, geometric distortion, and calibration in ab¬ 
solute flux. The absolute calibration factors that convert digital 
units into Wm~^nm~^sr~^ (referred to the central wavelength of 
each Alter) were recently revised in the OSIRIS pipeline and are 
slightly updated compared with those used before the Rosetta 
hibernation; they were successfully tested on Vesta (Eornasier et 



Phase Angle [deg] 

Fig. 1. Phase curve of the nucleus of 67P from OSIRIS observations in 
the NAC green filter (F23). The continuous line represents the HG lAU 
model that best fits the data, and the linear fit to the data is plotted as 
the dashed-dotted line. 


al., 2011). The absolute calibration factors were derived from 
Vega observations acquired during the instrument calibration 
campaign on 18 May 2014, and computed using the Vega and 
the Sun flux standard spectra from the HST CALSPEC catalog 
(www.caha.es/pedraz/SS S/HS TJOALS PEC). 

Eor each image, we computed the radiance factor (also 
known as //F) for each pixel. 


Radiance Factor(A) = 


7il(i, e, a, A) 


( 1 ) 


where I is the observed scattered radiance, the incoming so¬ 
lar irradiance at the heliocentric distance of the comet, and /, 
e, and a are the incidence, emission, and phase angles, respec¬ 
tively. The solar irradiance Fj_ is wavelength dependent and was 
calculated at the central wavelength of each Alter to be consistent 
with the methodology applied to derive the absolute calibration 
factors. 


3. Photometric properties from disk-averaged or 
integrated photometry 

3.1. Phase function 

Eor an irregular body such as the nucleus of comet 67P, the defi¬ 
nition of the disk-integrated phase function becomes ambiguous 
(Li et al., 2004) because the observed disk-integrated reflectivity 
depends upon the solar phase angle a and the illuminated sur¬ 
face seen during the observation. To built the phase function, we 
used the OSIRIS observations obtained with the NAC green Al¬ 
ter, centered at 535 nm. Rosetta approached the comet very fast 
in July-August 2014, with a resulting rapid increase of the spa¬ 
tial resolution (see Table [T]). The classical aperture photometry 
with a constant radius can therefore not be applied to evaluate 
the comet flux. We decided to integrate the comet signal over the 
pixels with a flux higher than 0.2x < F >, where < F > is the 
mean flux of the illuminated surface of the nucleus in a given 
image that was previously cut around the comet to avoid any po¬ 
tential background contributions (stars or cosmic rays). In this 
way, we integrated the comet flux over an area that is slightly 
larger than the projected surface of the comet obtained from the 
shape model for a given observing sequence, thus including the 
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Table 1. Observing conditions for the OSIRIS data. The filter names are followed by their central wavelength in nm at the first time that they 
appear in the table. The time refers to the start time of the first image of each sequence. A is the distance between the Rosetta spacecraft and comet 
67P. 


Camera 

Time 

filters 

phase (°) 

A (km) 

res. (m/px) 

NAC 

2014-07-21T14.34.16 

F22 (649.2), F23 (535.7), F24 (480.7), F16 (360.0), 

5.3 

4952 

93 

NAC 

2014-07-21T15.50.16 

F27 (701.2), F28 (743.7), F41 (882.1), F51 (805.3), 

F61 (931.9), F71 (989.3), F15 (269.3), & WAC-F71 (325) 

F22, F23, F24, F16, F27, F28, F41, F51, F61, F71, F15, & WAC-F71 

5.3 

4915 

93 

NAC 

2014-07-21T17.06.16 

F22, F23, F24, F16, F27, F28, F41, F51, F61, F71, F15, & WAC-F71 

5.2 

4878 

92 

NAC 

2014-07-21T18.22.16 

F22, F23, F24, F16, F27, F28, F41, F51, F61, F71, F15, & WAC-F71 

5.1 

4842 

91 

NAC 

2014-07-21T19.38.16 

F22, F23, F24, F16, F27, F28, F41, F51, F61, F71, F15, & WAC-F71 

5.0 

4805 

91 

NAC 

2014-07-21T20.54.16 

F22, F23, F24, F16, F27, F28, F41, F51, F61, F71, F15, & WAC-F71 

5.0 

4768 

90 

NAC 

2014-07-21T22.10.16 

F22, F23, F24, F16, F27, F28, F41, F51, F61, F71, F15, & WAC-F71 

4.9 

4731 

89 

NAC 

2014-07-21T23.26.16 

F22, F23, F24, F16, F27, F28, F41, F51, F61, F71, F15, & WAC-F71 

4.8 

4695 

88 

NAC 

2014-07-22T00.42.16 

F22, F23, F24, F16, F27, F28, F41, F51, F61, F71, F15, & WAC-F71 

4.8 

4658 

88 

NAC 

2014-07-25T08.11.21 

F22, F23, F24, F16, F28, F41, F71, & WAC-F71 

2.2 

3092 

58 

NAC 

2014-07-25T09.34.58 

F22, F23, F24, F16, F28, F41, F71, & WAC-F71 

2.2 

3074 

58 

NAC 

2014-07-25T11.00.58 

F22, F23, F24, F16, F28, F41, F71, & WAC-F71 

2.2 

3055 

58 

NAC 

2014-07-25T12.26.58 

F22, F23, F24, F16, F28, F41, F71, & WAC-F71 

2.2 

3037 

57 

NAC 

2014-07-25T13.52.58 

F22, F23, F24, F16, F28, F41, F71, & WAC-F71 

2.1 

3019 

57 

NAC 

2014-07-25T15.18.58 

F22, F23, F24, F16, F28, F41, F71, & WAC-F71 

2.1 

3000 

57 

NAC 

2014-07-25T18.10.58 

F22, F23, F24, F16, F28, F41, F71, & WAC-F71 

2.0 

2963 

56 

NAC 

2014-07-25T19.55.58 

F22, F23, F24, F16, F28, F41, F71, & WAC-F71 

2.0 

2940 

55 

NAC 

2014-07-28T22.10.25 

F27, F28, F41, F51, F61, F71, & WAC-F71 

1.3 

1984 

37 

WAC 

2014-07-28T22.10.25 

F18 (613), F71 (325), F13 (375), F14 (388), F31 (246), 

1.3 

1984 

195 

WAC 

NAC 

2014-07-28T22.10.25 
2014-07-28T23.25.31 

F41 (259), F51 (295), F61 (309), F81 (336) 

F22, F23, F24, F16, F27, F28, F41, F51, F61, F71, F15, & WAC-F71 

1.3 

1968 

37 

WAC 

2014-07-28T23.25.31 

F18, F71, F13, F14, F31, F41, F51, F61, F81 

1.3 

1968 

195 

NAC 

2014-07-29T00.45.31 

F22, F23, F24, F16, F27, F28, F41, F51, F61, F71, F15 

1.3 

1951 

37 

WAC 

2014-07-29T00.45.31 

F18, F71, F13, F14, F31, F41, F51, F61, F81 

1.3 

1968 

195 

NAC 

2014-08-01T11.50.14 

F82 (649.2), F23, F24, F27, F28, F41, F71, & WAC-F71 

9.0 

826 

16 

NAC 

2014-08-01T13.20.42 

F82, F23, F24, F27, F28, F41, F71, & WAC-F71 

9.3 

807 

15 

NAC 

2014-08-01T14.43.48 

F82, F23, F24, F27, F28, F41, F71, & WAC-F71 

9.7 

790 

15 

NAC 

2014-08-01T16.08.14 

F82, F23, F24, F27, F28, F41, F71, & WAC-F71 

10.1 

772 

14 

NAC 

2014-08-01T17.26.34 

F82, F23, F24, F27, F28, F41, F71, & WAC-F71 

10.4 

756 

14 

NAC 

2014-08-01T18.31.34 

F82, F23, F24, F27, F28, F41, F71, & WAC-F71 

10.8 

742 

14 

NAC 

2014-08-01T19.37.34 

F82, F23, F24, F27, F28, F41, F71, & WAC-F71 

11.1 

729 

14 

NAC 

2014-08-01T20.43.34 

F82, F23, F24, F27, F28, F41, F71, & WAC-F71 

11.5 

715 

13 

NAC 

2014-08-03T00.21.16 

F82, F23, F24, F27, F28, F41, F71, & WAC-F71 

28.0 

387 

7.3 

NAC 

2014-08-03T01.21.16 

F82, F23, F24, F27, F28, F41, F71, & WAC-F71 

29.1 

377 

7.1 

NAC 

2014-08-03T02.21.16 

F82, F23, F24, F27, F28, F41, F71, & WAC-F71 

30.2 

366 

6.9 

NAC 

2014-08-03T03.21.16 

F82, F23, F24, F27, F28, F41, F71, & WAC-F71 

31.6 

356 

6.7 

NAC 

2014-08-03T04.21.16 

F82, F23, F24, F27, F28, F41, F71, & WAC-F71 

33.0 

346 

6.5 

NAC 

2014-08-03T05.21.16 

F82, F23, F24, F27, F28, F41, F71, & WAC-F71 

34.4 

336 

6.3 

NAC 

2014-08-03T06.21.16 

F82, F23, F24, F27, F28, F41, F71, & WAC-F71 

35.9 

326 

6.1 

NAC 

2014-08-03T13.39.14 

F82, F23, F24, F27, F28, F41, F71, & WAC-F71 

40.7 

287 

5.4 

NAC 

2014-08-03T14.39.14 

F82, F23, F24, F27, F28, F41, F71, & WAC-F71 

40.7 

284 

5.3 

NAC 

2014-08-03T15.39.14 

F82, F23, F24, F27, F28, F41, F71, & WAC-F71 

40.6 

281 

5.3 

NAC 

2014-08-03T16.39.14 

F82, F23, F24, F27, F28, F41, F71, & WAC-F71 

40.6 

278 

5.2 

NAC 

2014-08-03T17.39.14 

F82, F23, F24, F27, F28, F41, F71, & WAC-F71 

40.6 

275 

5.2 

NAC 

2014-08-03T18.39.14 

F82, F23, F24, F27, F28, F41, F71, & WAC-F71 

40.6 

273 

5.1 

NAC 

2014-08-03T19.39.14 

F82, F23, F24, F27, F28, F41, F71, & WAC-F71 

40.6 

270 

5.1 

NAC 

2014-08-03T20.39.14 

F82, F23, F24, F27, F28, F41, F71, & WAC-F71 

40.6 

267 

5.0 

NAC 

2014-08-05T23.19.14 

F22, F23, F24, F27, F28, F41, F71, & WAC-F71 

48.9 

123 

2.3 

NAC 

2014-08-06T00.19.14 

F22, F23, F24, F27, F28, F41, F71, & WAC-F71 

49.5 

121 

2.3 

NAC 

2014-08-06T01.19.14 

F22, F23, F24, F27, F28, F41, F71, & WAC-F71 

50.0 

119 

2.2 

NAC 

2014-08-06T02.19.14 

F22, F23, F24, F27, F28, F41, F71, & WAC-F71 

50.8 

117 

2.2 

NAC 

2014-08-06T04.19.13 

F22, F23, F24, F27, F28, F41, F71, & WAC-F71 

52.2 

113 

2.1 

NAC 

2014-08-06T05.19.14 

F22, F23, F24, F27, F28, F41, F71, & WAC-F71 

53.0 

111 

2.1 

NAC 

2014-08-06T06.19.14 

F22, F23, F24, F27, F28, F41, F71, & WAC-F71 

53.9 

no 

2.1 


limb contribution as well. We evaluated the coma contribution to 
be 0.1-0.2 % of the comet signal in the green filter for the obser¬ 
vations up to 3 August, and 0.2-0.5 % for observations on 6 Au¬ 
gust, with the strongest coma contribution corresponding to the 
observation made at 04:19, when a jet departing from the Hapi 


region (see Thomas et al., 2015, for the region nomenclature for 
comet 67P) is clearly visible in the stretched image (see Fig.fTOl). 
The coma contribution is thus negligible given the low activity of 
the comet at these large heliocentric distances, and smaller than 
the uncertainties from absolute calibration, which are lower than 
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1% for the green filter (Tubiana et al., 2015a). Nevertheless, we 
assumed a conservative approach for flux uncertainties, consid¬ 
ering an error of 1.5% of the flux to take into account the coma 
contribution and the absolute calibration uncertainties. We then 
corrected the measured comet flux to take into account the fact 
that the comet spectrum is redder than that of the Sun, using 


£Fe(A)T(A)dA 
X Fco„,e>0)T(A)dA’ 


( 2 ) 


where and are the corrected and uncorrected cometary 
fluxes at the central wavelength Ac of the green filter, T(/l) is 
the system throughput (telescope optics and CCD quantum ef¬ 
ficiency), and Fo(/l) and Fcomet W are the solar (from the HST 
catalog) and cometary spectra (from Tubiana et al., 2011), re¬ 
spectively, both normalized to unity at the Ac of the green filter. 
The visual magnitude of the comet was reduced to the standard 
Bessel V filter by 


my = -2.5 



Fcomet{'F)Ty{A)dA \ + C, 


( 3 ) 


where Ty is the the transmission of the standard V filter taken 
from Bessel (1990). We determined the C constant as the dif¬ 
ference between the apparent solar magnitude of -26.75 (Cox, 
2000) and the convolution of the solar spectrum with the V fil¬ 
ter as in Eq.[3l The absolute magnitude V(l,l,Qf) reduced to the 
V-Bessel filter is shown in Fig.[T] 

We used the HG system (Bowell et al., 1989) to characterize 
the phase-function behavior of the nucleus and its absolute mag¬ 
nitude, Hy(l, 1,0). Our best-fit values for both parameters are G=- 
0.13+0.01, and F^y(l, 1,0)= 15.74+0.02 mag, which were deter¬ 
mined from chi-squared fitting. The nucleus displays a strong 
opposition effect: its linear slope (calculated for a > 1^) is p 
= 0.047+0.002 mag/°, and the magnitude from the linear slope 
(thereby excluding the opposition effect) is Hy//„ = 16.16+0.06 
mag. The scatter of the measurements obtained at similar phase 
angles (Fig.[T]) is related to the rotational phase of the comet. 

Attempts to determine the phase function of the nucleus 
of 67P were made using ground-based observations as well. 
Tubiana et al. (2011) sampled the phase-angle range 0.5°-10° 
and concluded that a linear approximation better represents their 
measurements than the lAU-adopted phase function. They deter¬ 
mined a very steep magnitude dependence on phase angle, with 
linear phase coefficients in the range 0.061-0.076 mag/°. Lamy 
et al. (2007) combined HST (Lamy et al., 2006) and ESO-NTT 
observations (Lowry et al., 2006) and obtained a slope parameter 
G=-0.45, which implies a very steep brightness dependence on 
phase angle and a possible strong opposition effect. This result 
might be affected by changes in the observing geometry that are 
due to the orbital motion of the comet. Lowry et al. (2012) deter¬ 
mined a slope parameter G=0.11 + 0.12 and, by a linear fitting, 
yS = 0.059 + 0.006 mag/°. 

The phase function determined from the OSIRIS observations 
has a lower linear slope than the values previously obtained from 
ground-based and HST observations, simply because it is evalu¬ 
ated for phase angles > 7^, excluding the contribution from the 
opposition effect. If we consider a linear fit of the OSIRIS data in 
the 1-10° range, then the linear slope is 0.074 mag/°, thus com¬ 
patible with previous measurements. Ciarniello et al. (2015) used 
data from the VIRTIS instrument and found a steeper slope for 
a < 15^ (0.082+0.016 mag/°) and a lower value (0.028+0.001 
mag/°) than we found for 25 < a < 110°. The discrepancies of 
the results between OSIRIS and VIRTIS data are related to the 



Wavelength [nm] 

Fig. 2. Average spectrophotometry of NAC and WAC observations ob¬ 
tained on 28 July, 23:25-23:45 UT. The continuous gray line represents 
a ground-based spectrum of 67P taken in 2007, when the comet was 
inactive (Tubiana et al., 2011). 


different phase-angle range used to compute the slope and to the 
fact that Ciarniello et al. (2015) derived the integrated magnitude 
from the disk-averaged albedo assuming that the nucleus of 67P 
is spherical with an equivalent radius of 1.72 km, while we di¬ 
rectly measured the integrated flux. 

The linear slope we determined is close to the average value 
for Jupiter-family comets (j8 = 0.053+0.016 mag/°, Snodgrass 
et al., 2011), very similar to the value found for comets Hartley 
2 (j3=0.046 mag/°, Li et al., 2013), Tempel 1 (j3=0.046 mag/°, 
Li et al., 2007a, 2013), and Borelly (j3=0.043 mag/°, Li et al., 
2007b), and similar to that of low-albedo asteroids (Belskaya & 
Shevchenko 1999). 

3.2. Disk-averaged spectrophotometry 

67P has been observed with several NAC filters since the very 
first resolved images, allowing us to investigate the spectropho- 
tometric properties of the nucleus. We present in this section the 
global spectrophotometric properties derived from the NAC and 
WAC images obtained on 28 and 29 July, 2014, when Rosetta 
was at a distance of about 2000 km from the comet, at the small¬ 
est phase-angle reached during the 2014 observations (Qf=1.3°). 
The corresponding spatial resolution was 37 m/px for the NAC 
(Table [B. The data set includes sequences that consist of ob¬ 
serving blocks of 11 NAC filters and 2 WAC filters repeated 
nine times during the rotational period of the comet (from 28 
July 14:39 UT to 29 July 00:45 UT ), but, as a result of pointing 
uncertainties, only the last three sequences, starting from 22:19 
UT, are centered on the comet and useful for our analysis. In 
addition, and only for this run, we obtained parallel observing 
sequences including all the 14 WAC filters repeated eight times 
from 28 July 15:35 UT to 29 July 00:55 UT. Unfortunately, the 
WAC images acquired with filters covering T > 500 nm are sat¬ 
urated or have fluxes beyond the linearity range of the detector, 
therefore we decided to discard them from our analysis. We de¬ 
rived a mean I/L value for each NAC and WAC image by inte¬ 
grating the cometary signal for all the pixels with a flux > 0.5x 
< I/F >, where < I/F > is the mean radiance factor of the illu¬ 
minated pixels on the cometary surface, and finally dividing by 
the projected surface of the comet derived from the stereopho- 
togrammetric (SPG) shape model (Preusker et al., 2015). 
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The spectrophotometry derived from the WAC and NAC images 
taken on 28 July, 23:25 UT, is shown in Fig.O with a superposed 
ground-based spectrum of 67P taken from Tubiana et al. (2011). 
The OSIRIS spectrophotometry in the visible and near-infrared 
range is fully consistent with ground-based observations, as also 
summarized in Tabled 

In general, the nucleus has similar spectrophotometric proper¬ 
ties in the NUV-VIS-NIR range to those of bare cometary nuclei 
(Lamy et al., 2004), of primitive D-type asteroids such as Jupiter 
Trojans (Fornasier et al., 2007, 2004), and of the moderately red 
Transneptunians and Centaurs (Fornasier et al., 2009). No clear 
absorption bands are visible at the spectral resolution of the fil¬ 
ters we used, and the spectral slope may change beyond 750 nm, 
being fainter in the NIR range than in the VIS. However, this 
change of slope might not be real because the 700 nm and 743 
nm fiuxes are contaminated by the coma emissions present at 
this large heliocentric distance (3.65 AU), as we show in Sect. 
4.3. 

Figure [2] shows that the fiux rises up in the mid-UV (240- 
270 nm range) for the two WAC filters F31 and F41 and for 
the NAC-F15. These two WAC filters, designed to study the CS 
emission and the adjacent continuum at 245 nm, are very diffi¬ 
cult to calibrate because they are affected by pinholes (Tubiana 
et al., 2015a), therefore some additional light may contaminate 
the signal. However, we checked the background signal for the 
two filters, and it is not as high as expected for a significant fiux 
contribution from pinholes defects. If we assume that this signal 
excess is at least partially real, then a tentative explanation of 
this broad band, which has a minimum close to 290 nm, may be 
SO 2 frost. A similar feature was seen on Jupiter satellites, for in¬ 
stance, the trailing hemisphere of Europa, where it was attributed 
to sulphur dioxide ice (Noll et al., 1995). Another similar feature 
was observed on comet Bowell 1980b (AHearn et al., 1984) and 
on the Centaur 2060 Chiron; Brosch (1995) attributed it to water 
ice or other frost such as NH 3 . 

Laboratory experiments on cryogenic films of SO 2 and SO 2 /H 2 O 
show that the SO 2 absorption band, usually located close to 280 
nm, is irreversibly temperature dependent (Hodyss et al., 2013) 
and it moves toward longer wavelength with the increase of the 
temperature. Hodyss et al. (2013) found a minimum close to 287 
nm for T ~ 140 K, and this minimum must slightly shift to longer 
wavelengths for higher temperatures, such as those measured on 
67P by the VIRTIS instrument, that is, 180-230 K (Capaccioni 
et al., 2015). For T > 160 K, the SO 2 can sublimate and thus the 
gas sulphur dioxide will be photolyzed into SO-hO and S - 1 - O 2 
with relatively high photodissociation coefficients (in the range 
of 1 X 10“"^ and 5 x 10“^ s“^ at 1 AU). It must be noted that 
the ROSINA instrument has identified sulphur compounds in the 
coma of 67P (Altwegg et al., 2015), and that the 67P UV spec¬ 
tra acquired with ALICE are compatible with the presence of 
SO 2 ice (Feaga et al., 2015). Therefore, the absorption feature at 
~ 290 nm could be due to SO 2 or any other ice containing that 
molecule in its structure. 

We also see an excess of the fiux in the OH filter centered at 
309 nm compared to the nearby continua (Fig. 13. To understand 
whether this higher fiux at 309 nm is due to the resonance fiuo- 
rescence of OH produced by water photolysis, we have, on one 
hand, studied the radial profiles of the emission both at 309 nm 
and in the nearby continua, and on the other hand, studied the 
remaining OH emission after subtracting the underlying contin¬ 
uum measured from the observation in the UV295 filter. 

From the study of the radial profiles measured from the observa¬ 
tions at 295 nm and 309 nm, no conclusive results can be drawn 
because both profiles show a similar behavior, which prevents 


concluding that the behavior of OH is characteristic of a gas 
species produced by a parent species. We also produced pure OH 
gas maps by subtracting the underlying continuum contribution. 
The resulting OH image is such that the emission excess is on 
the same order as the background signal. However, although the 
analysis is not conclusive because the images obtained on 28-29 
July were not designed to study gas emissions, it is worth noting 
that H 2 O in the gas phase was detected since early June 2014 
by the MIRO instrument (Gulkis et al., 2015), showing periodic 
variations related to nucleus rotation and shape. OSIRIS already 
detected an outburst of activity in April 2014 (Tubiana et al., 
2015b), and the ROSINA instrument (Altwegg et al., 2015) also 
measured H 2 O. This implies that the OH radical, coming from 
water photodissociation, is probably also present in the coma. 
Figure [2| also shows a slight increase of the I/F value at wave¬ 
lengths where CN fiuorescent emission takes place (i.e., 387 
nm^ A similar analysis to the one made for the OH case indi¬ 
cates that, if any CN surrounds the nucleus, the emission is on 
the same order of magnitude as that of the background. Thus, a 
clear confirmation of CN presence in the coma of 67P at these 
heliocentric distances cannot be drawn from our data. 

3.3. Phase reddening 



Phase angle [deg.] 

Fig. 3. Variation of the mean spectral slope of the nucleus of 67P, eval¬ 
uated in the 882-535 nm range and expressed in %/(100 nm), with the 
phase angle. 

We computed the disk-averaged spectral slope (S) following 
the definition of Delsanti et al. (2001) as 

^ _ _ ^882 - R535 _ 

^535 X (882 nm - 535 nm), 

where R 882 , and R 535 are the mean comet I/F values in the fil¬ 
ters centered on 882 nm and 535 nm. The mean spectral slope 
evaluated in the 882-535 nm range is shown in Fig [3] The ob¬ 
tained spectral slope varies between 11%/(100 nm) at a phase 
angle of 1.3° and 16%/(100 nm) at a phase angle of 52°, im¬ 
plying a significant phase reddening. This is the first time that 
significant phase-reddening effects are detected on a comet nu¬ 
cleus. Strong phase-reddening effects on comet 67P are also seen 
from the VIRTIS data both in the visible and near-infrared range 
(Ciarniello et al., 2015). Li et al. (2007a) used a different spectral 
slope definition, that is, it is normalized in the center of the wave¬ 
length range, and found a red slope of 12.5+1.0 % for comet 
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Table 2. Spectral slope of the nucleus of 67P, according to the spectral slope definition of Delsanti et al. (2001), derived from the literature and the 
data presented here. 


Date 

Slope (%/lOOnm) 

Wavel. range (nm) 

an 

Reference 

April 2004 

12 + 2 

545 - 797 

2.1 

Tubiana et al., 2008 

June 2004 

5 + 3 

545 - 797 

10.3 

Tubiana et al., 2008 

May 2006 

11+2 

545 - 797 

0.5-1.3 

Tubiana et al., 2008 

May 2006 

11+2 

500 - 850 

0.5-L3 

Tubiana et al., 2008 

July 2007 

10 + 2 

545 - 797 

6.0 

Tubiana et al., 2011 

July 2007 

11 + 1 

500 - 850 

6.0 

Tubiana et al., 2011 

July 2007 

10+1 

550 - 650 

6.0 

Lowry et al., 2012 

21 July 2014 

11.2+0.3 

535-882 

4.8-5.4 

This work 

25 July 2014 

11.3+0.3 

535-882 

2.2 

This work 

28-29 July 2014 

11.5+0.5 

535-882 

1.3 

This work 

1 Aug. 2014 

12.2+0.3 

535-882 

9.0-11.5 

This work 

3 Aug. 2014 

14.7+0.7 

535-882 

28-36 

This work 

3 Aug. 2014 

15.6+0.9 

535-882 

40-41 

This work 

6 Aug. 2014 

15.8+0.4 

535-882 

50-54 

This work 


Tempel 1 from the Deep Impact mission observations at a phase 
angle of 63° and no sign of phase reddening in the two high- 
phase angles they sampled (63° and 117°). With the same spec¬ 
tral slope definition, Li et al. (2013) measured a spectral slope 
of 7.6+3.6%/(100 nm), and of 8.6+10.6%/(100 nm) for comet 
Hartley 2 at phase angles of 85° and 93°, respectively. As a result 
of the large uncertainties in the spectral slope values, no conclu¬ 
sions on the phase reddening of Hartley 2 can be drawn. 

If we compute the spectral slope applying the same method as 
Li et al. (2007, 2013), that is, normalizing the reflectance value 
at the central wavelength range instead that at 535 nm, we And 
a slope value of 12.5%/(100 nm) at 52° that is similar to that of 
comet Tempel 1. 

The linear flt that best matches our data has a slope y = 
0.104+0.003 X 10“"^ nm“V^, and the disk-averaged spectral 
slope is estimated to be 11.3+0.2 %/(100 nm) at zero phase an¬ 
gle. 

The phase reddening was first noted by Gehrels et al. (1964) on 
the lunar surface, and it is potentially attributed to the increased 
contribution of the multiple scattering at large phase angles as 
the wavelength and albedo increase, plus a contribution of sur¬ 
face roughness effects (Sanchez et al., 2012; Hapke et al., 2012; 
Schroder et al., 2014). As reported in Fig.O the mean spectral 
slope varies not only with the phase angle, but, at similar phase 
angle, also with the rotational phase. Higher spectral slope val¬ 
ues correspond to redder surfaces, that is, when bluer regions like 
Hapi are not visible. For example, the peak in the spectral slope 
at Of=32-33° corresponds to color sequences taken on 3 August 
at UT 03:21-04:22 (Fig.O when the Hapi region, which is the 
bluest region of the comet (see Sect. 4.2), was not visible at all. 

3.4. Hapke model on disk-averaged photometry 


study the photometric properties of the comet’s surface. We used 
the Hapke model (Hapke, 1993), a semi-theoretical assembly 
of functions that describes the contribution of single-scattering, 
multiple-scattering, opposition effect and roughness to repro¬ 
duce the measured reflectance. The model has experienced sev¬ 
eral improvements with more physical description since its in¬ 
troduction in the 1980s (Hapke, 1993, 2002, 2008, 2012), which 
is the reason that there are several different versions. For a comet 
that is expected to be composed of dark phases rich in carbon, the 
coherent-backscattering mechanism (Shevchenko & Belskaya, 
2010; Shevchenko et al., 2012) is not expected to play a ma¬ 
jor role. For the single-particle-phase function, we adopted the 
single-term of the Henyey-Greenstein (HG) form (Hapke 2012), 
with one asymmetry parameter, g. A negative g value represents 
backscattering, a positive one forward scattering, and g = 0 
stands for isotropic scattering. 

Since the first resolved images, the nucleus of 67P presented 
a highly irregular shape. However, Li et al. (2003) have re¬ 
ported that the influence of shape is irrelevant when obtaining 
the Hapke parameters when the data are flt with phase angles 
lower than 60°, which is our regime. Therefore, we adopted 
a five-parameter version of the model, built for spherical bod¬ 
ies, that includes the single-scattering albedo the asymmetry 
factor g^, the average roughness angle 0, and the amplitude 
and width hg of the shadow-hiding opposition surge. The disk- 
integrated Hapke expression for the radiance factor can be writ¬ 
ten as 

^ = ^(a.^)[( Y [d + BsH,(a))PHg(a,g,) - l] + y (1 - ro,)) 

^1 - sin ^ tan ^ In |cot j + a + (t: - a) cos Qf)|, (4) 


We computed the disk-averaged reflectance in different Al¬ 
ters by integrating the cometary signal for all the pixels with a 
flux > 0.5X < I/F >, where < I/F > is the mean radiance fac¬ 
tor of the illuminated pixels on the cometary surface, and Anally 
dividing by the projected surface of the comet derived from the 
SPG shape model. 

The disk-averaged multiwavelength phase curve (in radiance 
factor, I/F) in the L3°-54° phase-angle range allows us to 


where a is the phase angle, A is the wavelength were the j 
is measured, and ro is defined as 


1 - yi - Wj 

1 + yi - ’ 


( 5 ) 
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Table 3. Hapke parameters derived from disk-averaged reflectance for different Alters. 


/l(nm) 

Wa 

9a 

Bo 

hs 

6 

A 

^geo 

^bond 

325 

0.028 + 0.001 

-0.35 + 0.03 

1.83 + 0.03 

0.029 + 0.005 

15 

0.0316 + 0.0030 

0.0088 + 0.0007 

480 

0.035 + 0.001 

-0.43 + 0.03 

1.91 + 0.09 

0.021 + 0.004 

15 

0.0554 + 0.0024 

0.0119 + 0.0001 

535 

0.037 + 0.002 

-0.42 + 0.03 

1.95 + 0.12 

0.023 + 0.004 

15 

0.0589 + 0.0034 

0.0123 + 0.0001 

649 

0.045 + 0.001 

-0.41 + 0.03 

1.97 + 0.09 

0.026 + 0.007 

15 

0.0677 + 0.0039 

0.0157 + 0.0001 

700 

0.050 + 0.001 

-0.38 + 0.02 

2.22 + 0.06 

0.027 + 0.002 

15 

0.0720 + 0.0031 

0.0173 + 0.0001 

743 

0.053 + 0.004 

-0.40 + 0.03 

1.96 + 0.07 

0.023 + 0.003 

15 

0.0766 + 0.0031 

0.0178 + 0.0001 

882 

0.052 + 0.003 

-0.40 + 0.04 

2.08 + 0.09 

0.027 + 0.006 

15 

0.0780 + 0.0038 

0.0179 + 0.0001 

989 

0.066 + 0.005 

-0.35 + 0.03 

2.07 + 0.05 

0.032 + 0.004 

15 

0.0820 + 0.0039 

0.0223 + 0.0001 


where w is the single scattering albedo. The shadow-hiding op¬ 
position effect BsH^(a) and the shadowing function K;i(a, 6) are 
given by Hapke (1993): 


Bo, hs) 


Bo 


1 -f 


tan h 

~h7 


( 6 ) 


f 

a 

l -0329 

tan 6 tan — 

[ 2J 


1/2 _ _ Of 

- 0.52^tan ^ tan — 
2 


(7) 


and the single-term Henyey-Greenstein function is expressed as 




(1 + 20-^008 O'+ g2)3/2- 


( 8 ) 


To model the observed data with the Hapke formula, we de¬ 
veloped a procedure that consists of randomly picking 30 initial 
conditions for each parameter for a given filter, and we searched 
for the global minimum between the observations and the 
fit using Broyden-Fletcher-Goldfarb-Shanno (BFGS) algorithm 
(Byrd et al., 1995) in the Basin-Hopping method (Wales & Doye, 
1997). The minimum was searched for inside the boundaries 
wa = {0.02,0.1}, qa = 1-1.0,1.01, ^0 = (0.2,2.5}, and/z, = 
(0.0,0.5}. Preliminary tests showed that 6 is ill constrained due 
to the lack of points at phase angles larger than 60°, so we de¬ 
cided to fix it to 15°, the value determined from the model of 
Hapke (2002) on the disk-resolved images, as discussed in Sect. 
4.2. 

The Basin-Hopping method has been used to solve com¬ 
plex molecular systems (e.g., Prentiss et al., 2008; Kim et al., 
2014) and is based on scanning several local minima to esti¬ 
mate the global minimum. When the solution was found, we 
ran a Levenberg-Marquardt algorithm of the MINPACK library 
(More, 1978) for each parameter individually to adjust them to a 
better precision. The best solution was chosen from the 30 tests 
as the test that returned the lowest value. The standard de¬ 
viation of the parameters was calculated from the spread of the 
solutions of each test. Because we used a method that estimates 
the global solution, deviations among the solutions found give 
an estimate of the error bars of the different parameters. Under¬ 
sampling and correlations among the Hapke parameters (Bowell 
et al., 1989) are the main sources of errors when flux uncertain¬ 
ties are negligible. 

The obtained phase curve shows a steeper behavior starting from 
~ a =5°, probably because of the shadow-hiding opposition 
effect. The obtained amplitude of the opposition effect is no 
smaller than 1.78, and its width is no lower than 0.023. The 
comet has thus a sharper opposition effect than that found from 



Wavelength [nm] 

Fig. 4. Geometric albedo for different Alters derived from Hapke mod¬ 
eling of the disk-averaged photometry. 


an average phase function for C-type asteroids (^o = 1-03 and 
hs = 0.025, Helfenstein & Veverka, 1989). 

The Hapke parameters we obtained for the 67P disk-average 
reflectance (Table O are very similar to those found by Li et al. 
(2012) on disk-averaged reflectance at 650 nm of comet Tern- 
pel 1 from the Deep Impact data (i.e., wesonm = 0.043 + 0.006, 
qa = -0.48 + 0.02, g = 20 + 5 and = 0.059 + 0.046 ). We 
report all multifllter Hapke solutions from the disk-averaged re¬ 
flectance in Table [3l Hapke modeling also enables deriving the 
geometric albedo, which is reported in Fig. |4] for the different 
wavelengths we investigated. The Hapke parameters show no 
clear wavelength dependence of qa, Bo and hs, within the un¬ 
certainties. This behavior is expected when the shadow-hiding 
effect is the main cause of the opposition surge, while for the 
asymmetry factor, it means that the average scattering path of 
light inside the surface particles is small and that light is rapidly 
absorbed. 

4. Disk-resolved photometry, colors, and 
spectrophotometry 

Along with the OSIRIS NAG images, we derived additional 
files that contain geometric information about the illumination 
and observation angles. For each NAG image pixel, we com¬ 
puted solar incidence, emission, and phase-angle information. 
This computation considers all relevant geometric parameters, 
such as the camera distortion model, the alignment of the cam¬ 
era to the Rosetta spacecraft, the orientation of the spacecraft 
(with reconstructed orbit position and pointing) with respect to 
the 67P nucleus, and a 3D shape model of the comet. The shape 
model and the improved spacecraft orientation data are based 
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Fig. 5. RGB images, in false colors, of the nucleus of 67P from different observing runs. The color images are produced using the filters centered 
on 480 nm, 649 nm, and 882 nm. The spatial resolution varies from 90 m/px on July 21 (at the top) to 2.1 m/px on 6 August (at the bottom). 
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2014-08-21 

15h40 


2014-08-22 

01h40 



Fig. 6. RGB images, in false colors, of the nucleus of 67P from images taken on 21 and 22 August with the filters centered on 480 nm, 649 nm, 
and 882 nm. Rosetta was at a distance of 70 km from the nucleus, and the corresponding resolution is 1.3 m/px. 


upon a stereophotogrammetric analysis (SPG) of OSIRIS NAC 
images. SPG has already been successfully applied to OSIRIS 
NAC images taken during the Rosetta flyby of the asteroid (21) 
Lutetia (Preusker et al., 2012) and to imaging observations of 
many other planetary missions. The SPG-based shape model of 
the nucleus of 67P used in this analysis consists of about 6 mil¬ 
lion vertices (about 12 million facets) and comprises a sampling 
distance of about 2 m on the nucleus surface (Preusker et al., 
2015). 

Color cubes of the surface were produced by stacking registered 
and illumination-corrected ((//F’)corr) images. For the illumina¬ 
tion correction, we used a simple Lommel-Seeliger disk law. 


D{i, e) = 


2 cos(0 

cos(0 -f cos(^) ’ 


(9) 


The Lommel-Seeliger law comes from radiative transfer theory 
when considering a single-scattering particulate surface (Fair- 
bairn, 2005). This law has been proved to be suitable for smooth 
low-albedo surfaces, such as C-type asteroids and cometary nu¬ 
clei, due to the predominance of single scattering, especially 
at low phase angle, as predicted by radiative transfer models 
(Hapke, 1981). Disk-corrected images are hence obtained as 


(IIF\orr(a,A) 


7Tl(i, e, a, 
D{U e) 


( 10 ) 


The nucleus of comet 67P has a complex shape, with differ¬ 
ent types of terrains and morphological features (Thomas et al., 
2015), and it shows some color variations, in particular in the 
Hapi region (Figs. [5] and O. This region is located between the 
two lobes of the comet, it is both the brightest and the most ac¬ 
tive surface on the comet at the large heliocentric distances mon¬ 
itored during the approach phase and bound orbits of Rosetta. 
Hapi is also the source of the most spectacular jets seen at he¬ 
liocentric distances > 3 AU (Lara et al., 2015). Other bright re¬ 
gions appear at the boundary between the Aten and the Ash area 
together with bright, circular features on the Seth unit. At higher 
resolution, several bright spots associated with exposed water ice 
are seen in all types of geomorphologic regions on the comet, 
and these features are fully discussed in Pommerol et al. (2015). 
The peculiar colors of some regions of the comet, especially for 
the Hapi region, are confirmed by the bluer spectral slope com¬ 
pared to the darker regions, as we discuss in Sect. 4.2. Moreover, 
the Hapi region itself shows difference in colors when observed 
at higher resolution (Fig. O, and these differences are observ¬ 
able through different Alters and within multiple phase angles 
and resolutions, hence they are independent of observational ge¬ 
ometries and shape corrections. 


4.2. Albedo map and disk-resolved photometry 


4.1. Color variations 

OSIRIS data show the 67P surface with unprecedented resolu¬ 
tion compared to other cometary nuclei visited by space mis¬ 
sions. Figures [5] and [6] show the RGB maps, in false colors, 
generated using NAC images acquired from 25 July to 22 Au¬ 
gust with the Alters centered on 882 nm, 649 nm, and 480 nm. 
Images were first co-registered, then RGB images were gener¬ 
ated using the STIFF code (Berlin, 2012), a software developed 
to produce color images from astronomic FITS data. 


We investigated the local photometry and reflectance variations 
of the observed surface of the nucleus at 649 nm using Hapke’s 
(2002) reflectance model. The disk-resolved data offer many 
combinations of emission, incidence, and azimuth angles, allow¬ 
ing us to constrain some parameters in a better way than using 
disk-integrated data. We used topographic information based on 
the SPG shape model, together with the observed reflectance at 
649 nm on all the illuminated pixels of the selected images de¬ 
scribed in Table [U In the Hapke (2002) model, the bidirectional 
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Fig. 7. Goodness plot of the modeled I/F with respect to measured I/F 
(at 649 nm) of Hapke (2002) modeling of the nucleus of 67R 



Normal Albedo at S49 nm 



DlOC QlO? OlOS 0l09 


Fig. 8. Normal albedo polar map of the nucleus of 67P in the orange 
filter (649 nm). The longitude 0^ in the Cheops reference frame is to the 
bottom (Preusker et al., 2015). 

reflectance (Ref) at a given incidence (i), emission (e), and phase 
angle (a) is described as 

CO Un 

Ref(i, e, a) = - - [Phg(a, g)BsH(g) + M(fio,fi)]BcB(g), 


where po = cos(i), p = cos(^), Phg is the single-term 

Henyey-Greenstein function, Bsh the shadow-hiding opposition 


term, Bcb the coherent-backscattering term, and M the multiple¬ 
scattering term described in Hapke (2002, Eq. 17). To provide 
equivalent weight to the observations obtained at different phase 
angles, all the images were first binned at an interval of about 2.0 
degrees into a 3D grid composed of the three geometric angles. 
The I/F pixels were averaged over each bin. In this way, we re¬ 
moved the influence of intrinsic albedo variations on the fitting. 
To avoid unfavorable observation geometries due to large inci¬ 
dence and emission angles, we Altered the data and considered 
pixels with incidence and emission angles lower than 70^. This 
constrain excludes extreme geometries near the limb or near the 
terminator. A least- ht using the Levenberg-Marquardt algo¬ 
rithm was performed to And the global Hapke parameters, giving 
an RMS error of 6%. The modeled I/F with respect to the mea¬ 
sured I/F values are shown in Fig. [71 The Hapke parameters that 
best fit our data are given in Tabled] 

To test the robustness of the results, we also com¬ 
puted the Hapke parameters using the latest version of the 
Hapke model (Hapke, 2012), which also takes into ac¬ 
count the top layer porosity, given by the porosity factor 
K = -ln(l - 1.2090^/^)7(1.2090^/^), where 0 is the Ailing 
factor, and the porosity is defined as 1-0. The complete equation 
for this version of the model can be found in Helfenstein & 
Shepard (2011). We At all images together using the Broyden- 
Fletcher-Goldfarb-Shanno (BFGS) algorithm (Byrd et al., 
1995), which solves nonlinear problems by approximating the 
first and second derivatives in an iterative procedure to search 
for a local minimum. 

The boundaries used for the parameters are = 

{0.02,0.1}, = 1-1.0,1.01, Bo = (0.2,2.5}, = (0.0,0.5}, 

and 6 = {5^,90^}. We decided against using the basin-hopping 
method in this case because of constraints in the convergence 
time. Nevertheless, we used the disk-integrated results as initial 
conditions to the BFGS algorithm, expecting that the best 
solution might be situated around the same minimum. 

At the end of the procedure, we observed that the coherent- 
backscattering terms have remained invariant to the final 
solution, so we can rule out the coherent-backscattering process 
as a significant effect to the opposition regime at least for 
a > 1.3°. The best solutions from Hapke (2012) modeling are 
reported in Table |4j We observe that the regolith’s particles are 
backscattering, with a g parameter value close to the one found 
from the disk-integrated analysis. The single-scattering albedo 
value is lower than the one found from the disk-integrated 
analysis or from the Hapke (2002) model. This is due to the 
inclusion of the porosity term (Hapke 2008, 2012) that also 
acts on the opposition effect and multiple-scattering formulae. 
To compare the single-scattering albedo from the two models, 
the wx from the Hapke (2012) solution must therefore be 
multiplied by the porosity factor K, which has a value of 1.2 
for a porosity of 0.87. In doing so, we obtained a Hapke (2012) 
porosity-corrected w value of 0.041, very close to the w value 
found from the Hapke (2002) model. 

We And a porosity value reaching 87% for the 67P nucleus 
top surface layer, which is higher than the 40-60% value 
found by the KOSI experiments (Sears et al., 1999; Kochan 
et al., 1998) that simulated cometary top surfaces. Faboratory 
experiments also showed the formation of a thick dust mantle 
on cometary-analog materials that may be exposed to dust 
settling and pulverization, facilitating the formation of a very 
weak and porous mantle (Mohlmann, 1995; Sears et al., 1999). 
Moreover, recent works have evoked fractal aggregates (Bertini 
et al., 2007; Fevasseur-Regourd et al., 2007; Fasue et al., 2010) 
as the best analogs of cometary dust. Surfaces composed of 
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Table 4. Results from the Hapke (2002, and 2012) modeling from disk-resolved images taken with the NAC orange filter centered on 649 nm. 


Model 

w 


Rose 

hg 


Rocb 

hcB 

Porosity 

K 

Geom. alb. 

Hapke 2002 

0.042 

-0.37 

2.5 

0.079 

15 

0.188 

0.017 

_ 

_ 

0.064 

Hapke 2012 

0.034 

-0.42 

2.25 

0.061 

28 

_ 

_ 

0.87 

1.2 

0.067 


such material have a porosity ranging from 80% to 90%, in 
agreement with what we found from the Hapke model. 

The single-scattering albedo, roughness angle, and geomet¬ 
ric albedo found from Hapke (2002, 2012) modeling are very 
similar to the values found for comets 9P/Tempel 1 and 81 P/Wild 
2 (Li et al., 2009, 2012). The asymmetry factor found for comet 
67P is lower than the one derived for comets Tempel 1, Hart¬ 
ley 2, and Wild 2 (Li et al., 2007a, 2009, 2013), implying that 
the regolith on comet 67P is less backscattering than for these 
comets. Ciarniello et al. (2015) analyzed disk-resolved data of 
comet 67P from VIRTIS observations and found an asymme¬ 
try factor of -0.42, a roughness of 19°, and a single-scattering 
albedo of 0.052+0.013 at 550 nm. While the 6 and g values 
they determined are similar to the one we found, the single¬ 
scattering albedo is much higher than the values we determined 
both from the disk-averaged and disk-resolved analysis of the 
OSIRIS data at 535-649 nm. It must be noted that Ciarniello 
et al. (2015) used a simplified version of the Hapke model, ne¬ 
glecting the opposition-effect terms and the multiple-scattering 
functions (implying an overestimation of the single-scattering 
albedo), which they applied to VIRTIS data covering the 27- 
111° phase-angle range. The different formalism and wavelength 
range covered explain the discrepancy in the single-scattering 
albedo values. 

If we compare the Hapke parameters from disk-resolved and 
disk-integrated photometry (Tables [3] and ID), we see that the in¬ 
tegrated solution gives a similar w but much lower hg values than 
the disk-resolved solutions. This discrepancy is related to the dif¬ 
ferent formalisms applied to compute the parameters, and to the 
fact that the effects of the different parameters are entangled in 
the disk-integrated function. Nonetheless, without data under 1.3 
degrees of phase angle, it is difficult to precisely determine the 
amplitude of the opposition effect. Therefore, Hapke parameters 
obtained from the disk-resolved data should be preferred as they 
are computed by modeling the full brightness variations across 
the comet surfaces. It is noteworthy that despite the discrepan¬ 
cies between disk-integrated and disk-resolved results, we ob¬ 
tained very similar geometric albedo values at 649 nm from the 
different solutions (Tables [3] and ID). 

The albedo map shown in Fig. [ 8 ] was then generated using a 
rotational movie sequence obtained on 3 August in which most 
of the comet surface was observed several times at very different 
incidences and emission angles. We used the parameters derived 
from the Hapke (2002) model to correct the refiectance in the 
following manner: 


-^corrected — ^^ measured 


^I^Hapkei}ref ■) ^ref-> (^ref) 
iI^Hapkeijobs’i ^obsr> ^obs^ 


( 12 ) 


where ///^measured Stands for the observed refiectance at the 
given incidence, emission, and phase angles (iobs^^obs^O'obs)^ 


while iref, eref, and are/ are the geometric angles at which the 
refiectance IlFcorrected is supposed to be. Here we have chosen 
/^gy^=0°, eref=0°, and aref=0°, so the corrected I/F should be the 
normal albedo. 

The //F’corrected of each surface element (facet) was then com¬ 
puted and stored as soon as a facet was illuminated. For each 
facet we obtained multiple measurements that were finally aver¬ 
aged to compute the normal albedo. Approximately 50 images 
were used to generate the albedo map represented in Fig.O The 
same approach was used by Leyrat et al. (2010) to study the 
photometric properties of asteroid Steins from the OSIRIS data 
acquired during the Rosetta flyby. 

The comet nucleus is dark and the estimated geometric 
albedo from the Hapke parameters is 6 .5+0.2% at 649 nm. This 
value is slightly higher but still similar to that found for comets 
Wild 2 and Tempel 1 (see Table 1 in Li et al., 2013 for a summary 
of the photometric parameters of comet nuclei visited by space 
missions). The histogram of the albedo distribution derived from 
the albedo map shown in Fig. [ 8 ] was fit with a Gaussian function 
that peaks at 0.063 and has a standard deviation cr of 0.01. Thus 
the surface reaches albedo variations of ~ 30% within 2cr of the 
albedo distribution. 

Some very bright spots in Fig. 8 seem to be related to high- 
frequency topographic features that are only poorly described by 
the shape model (i.e., surrounding terrains of the Imhotep region 
and cliffs between Hapi and Hathor). Overall, Hapi is the bright¬ 
est region, also showing some local albedo variations, as already 
noted in the RGB images (Fig. ( 6 ]). The Imhotep area also shows 
a significant albedo increase over a wide surface. Both the Hapi 
and Imhotep regions have a smooth surface, suggesting that the 
bright albedo may be at least partially due to the texture of the 
surface. On the other side, the Seth geological area is the darkest 
of the northern hemisphere. 

4.3. Local spectrophotometry and composition 

We computed the spectral slope for each pixel in the 535-882 
nm range on co-registered and illumination corrected images as 

^ _ Rcssijcr) -Fc535(a) 

R^ 535 (rr) X (882 nm - 535 nm), 

where Rc 882 (q^), and Rc 535 (Qf) are the comet I/F values in a 
given pixel, corrected for the illumination conditions using the 
Lommel-Seeliger disk law, in the filters centered on 882 nm and 
535 nm. Images were normalized to the green filter, centered 
on 535 nm, to be consistent with most of the literature data on 
primitive solar system bodies whose spectral slopes are usually 
computed on spectra normalized in the V filter or at 550 nm. 
We decided to correct the images for the illumination conditions 
using the simple Lommel-Seeliger disk law and not using the 
Hapke formalism simply because the Hapke parameters from 
disk-resolved images are only available for the filter centered on 


Article number, page 11 ofm 








A&A proofs: manuscript no. 25901_xarv 


l/F(a=1.3°) [%] 


4.6 4.8 5.0 5.2 5.4 5.6 5.8 6.0 
2001.I 



o r.1 

0 50 100 150 200 


Spectral slope 


8 9 10 11 12 13 14 15 

2001.J 



n r.1 

0 50 100 150 200 



Wavelength [nm] Wavelength [nm] 


Fig. 9. 1/F for the green filter centered on 535 nm, slope map, absolute and relative reflectance for some selected regions on the nucleus of 67P for 
the 29 July observations taken at 00:43 UT and at cr = 1.3°. The slope is computed in the 535-882 nm range, after normalization at 535 nm, and it 
is in %/(100 nm). The I/F map and the reflectance are corrected for the illumination conditions using the Lommel-Seeliger disk law. 


649 nm, and not at several wavelengths. 

Figure [9] shows the I/F map, illumination corrected using the 
Lommel-Seeliger disk law, of the comet obtained on 29 July 
with the NAC green filter, at the lowest phase angle (1.3°) sam¬ 
pled during the 2014 observations, as well as the spectral slope 
map, evaluated in the 882-535 nm range, and the relative and 
absolute reflectance of three regions selected over the surface. 
The disk-corrected image at 535 nm has a mean I/F value of 
5.2% at cr = 1.3°, and shows albedo variations up to -fl6% in 
the brightest Hapi region, and -9% in the darker Ash region. The 
reflectance is anticorrelated with the spectral slope, and in par¬ 
ticular, the bright Hapi region is also characterized by a bluer 
spectral slope than the darker medium-reflectance regions. We 
computed the absolute and relative spectrophotometry in a box 
3 px X 3 px wide (i.e., llOmx 110 m) on three selected regions 
representing dark, medium, and high reflectance areas (red star, 
blue squares, and black circles in Fig.O. The selected dark area 


shows a linear red slope without obvious spectral features, and it 
has a steeper spectral slope. The bright area selected in the Hapi 
region has a different spectral behavior, bluer than the other two 
regions, with a clear change of slope beyond 650 nm and possi¬ 
bly an absorption feature in the 800-900 nm range. All the three 
regions show an enhancement of the fluxes at 700-743 nm due 
to cometary emissions that we describe in the following para¬ 
graphs. 

We then studied the spectrophotometry of the comet obtained 
at higher spatial resolution during the observations on 6 August. 
Figure [TT] shows the relative reflectance, normalized at 535 nm 
and averaged over a box 5x5 pixels wide (~ 10 x 10 m^) for 
ten selected regions on the nucleus as reported in Fig. [TOl To 
study the relative spectrophotometry, we sampled three differ¬ 
ent groups of terrains, from the brightest and spectrally less red 
Hapi region to the reddest surfaces on the comet (Figs. [T0| and 
E]). The brightest regions, in particular Hapi, are clearly less red 
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Fig. 10. Three images acquired on 6 August at 00:19, 04:19, 06:19 UT, and again 04:19 UT, from left to right. The first three images from the 
left side show the selected regions where the spectrophotometry on the nucleus has been computed, as represented in Fig. nn The last image on 
the right side is taken at 04:19 UT but stretched to show the coma features. The points in this image represent different nucleus and coma regions 
whose flux is represented in Fig. cn 
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Fig. 11. Relative reflectance from ten selected regions as displayd in 
Fig. [T0| The dashed-dotted and dashed lines represent the mean spec¬ 
trophotometry for the TNO taxonomic classes BR and IR (Barucci et 
al., 2005), respectively. 


Fig. 12. Flux of the comet nucleus (divided by 100) and of the coma 
at several projected distances (the selected regions are represented in 
Fig- uni last image on the right) from the nucleus surface for the obser¬ 
vations on 6 August acquired at 04:19 UT. 


in term of spectral slope. We interpret this spectral behavior as 
being caused by a higher abundance of water ice at the surface, 
on the basis of the correlation between bluer slope and ice abun¬ 
dance found for comets Hartley 2 and Tempel 1 (Li et al., 2013; 
Sunshine et al., 2006, 2013), and on the basis of the results from 
experiments on ice or dust mixture sublimation (Pommerol et 
al., 2015). Moreover, the Hapi region is often shadowed by the 
Hathor region, it is not illuminated at perihelion due to seasonal 
effects, and thus it is a favored region to retain ices on its sur¬ 
face. It must be noted that no clear water ice absorption bands 
have been detected with the VIRTIS infrared-imaging spectrom¬ 
eter during the August 2014 observations at resolutions of 15- 
25 m/px (Capaccioni et al., 2015), even though the shape of the 
strong absorption beyond 3 yum in the Hapi region seems to be 
consistent with the presence of low amounts of water ice. The 
absence of large areas of water-ice-rich mixtures in VIRTIS data 
is attributed to the lower spatial resolution of these data than that 
of the OSIRIS images, together with the presence of dark-phased 
and non-volatile materials that may mask the water ice spectral 
absorptions. 

The reddest spectral regions are found in both lobes of the comet, 
in particular nearby the Hatmehit depression in the small lobe 
and in the Apis region in the large lobe. The different spectral be¬ 


haviors of the selected regions over the nucleus are comparable 
with those of other cometary nuclei, with Jupiter Trojan asteroids 
and with part of the transneptunian (TNO) objects (Fornasier et 
al., 2007, 2009; Lamy et al., 2004). We also plot in Fig. [TT] the 
two TNO taxonomic classes BR and IR (defined in Barucci et 
al., 2005) whose spectrophotometry is close to that observed in 
different areas of the comet. These two classes are intermediate 
between very blue and very red TNOs. 

With the spectral resolution given by the filters used it is not 
clear whether some absorptions features are present in the near- 
infrared regions, and the interpretation is complicated by the fact 
that we clearly see an excess of flux in the 700-750 nm region. 
To ascertain the origin of this flux excess, we studied the behav¬ 
ior of this emission at several locations both on the nucleus and 
in the coma. More concisely, for the images obtained on 6 Au¬ 
gust, we selected several locations on the nucleus, immediately 
(as projected on the sky) above the surface, at the inner limb, 
and in the coma to average the radiance in squares of 10 x 10 
m^, as was done for the analysis of different regions over the nu¬ 
cleus, for every filter. These measurements were made both in 
cometary regions where jet-like features are clearly visible and 
where they seem to be absent. 

Figure [T2| shows the coma radiance on the nucleus and at differ- 
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ent projected distances, from oservations on 6 August acquired 
at UT 04:19. Whereas at A < 649.2 nm, the I values on the nu¬ 
cleus and in the the dust coma follow the same trend, at 701 nm 
there is a noticeable increase in the inner coma locations, and the 
peak values decrease with increasing projected distance from the 
nucleus surface. A similar behavior is observed at 743 nm and 
tentatively ai A > 880nm, although the latter cannot be clearly 
confirmed as there are no data acquired with adjacent filters to 
allow any emission band isolation. On the directly illuminated 
nucleus surface both the band-like structure at 701-743 nm and 
the I/F strong increase in the near-IR cannot be detected when 
compared with the nucleus spectrum itself. That this behavior 
is mainly observed at coma positions indicates that this fiux is 
unrelated to the nucleus surface composition, but is related to 
the coma. The enhancement of the fiux in the 700-750 nm re¬ 
gion must then be associated with the cometary emissions in the 
coma. Furthermore, it was also verified that it can be clearly seen 
in shadowed areas, where the strong nucleus signal is absent. 
Using Mie-scattering models, we have ruled out that dust grains 
were responsible for the high values of I at 701 nm. Thus, the 
other possible candidates are emissions produced by gas species 
in the coma of 67P. From high-resolution spectra of comets 
(Cochran & Cochran, 2002), we conclude that the most likely 
candidates to give rise to such an enhanced emission are H 20 ^ 
in its (0,6,0) band at 700 nm, the (0,2,0)-(0,0,0) band at 743-745 
nm, and the band system (0,6,0) of NH 2 that overlaps with the 
water bands. Water, the parent species of H20^, was reported by 
Hassig et al. (2015), and Gulkis et al. (2015) from data acquired 
by the ROSINA and MIRO instruments, respectively, as early as 
August 4, 2014. The processes by which H 2 O is ionized once is 
still under study, although electronic impact and photoionization 
are among the most plausible ones. 

In addition to this coma feature, the spectrally bluer to moder¬ 
ately red regions show a rise of the fiux at 989 nm (Fig.[TT]), and 
this behavior was also seen in the VIRTIS spectra (Capaccioni 
et al., 2015). In particular, the spectrally bluer regions like Hapi 
seem to present a faint absorption feature centered between 800 
nm and 900 nm, which may be associated to the Fe^^ - Fe^^ 
charge transfer absorption band in silicates. Unfortunately, the 
low spectral resolution of the OSIRIS observations and the con¬ 
tamination of the adjacent filters in the 700-750 nm regions by 
coma emissions complicates the interpretation of this potential 
absorption feature. 

4.4. Correlations between albedo, color variations, and 
geology 

According to the albedo and color variations described in the 
previous paragraphs, the surface of the nucleus shows some 
variability. To investigate the correlations between geology and 
spectral properties, we computed the spectral slope maps (in 
%/100 nm, evaluated in the 882-535 nm range using the same 
method as described in Sect. 4.3) for the observations obtained 
on 1 and 6 August at different comet rotational phases. Fig¬ 
ures 13 and 14 show these spectral slope images, acquired when 
Rosetta was at a distance of about 800 km and 120 km from 
the comet, which yielded a spatial resolution of 14 m/px and 2 
m/px, respectively. The spectral slope was evaluated as the mean 
value over a 3x3 pixel squared box. These maps confirm the pe¬ 
culiarity of the Hapi region at increasing resolution scales. They 
also show the phase-reddening effect previously seen in the disk- 
averaged spectrophotometry, as the spectral slope of the same 
areas increases in the 6 August maps, when the phase angle was 
~ 50°, compared to the 1 August images, where the phase angle 


was 10°. 

Figures [13] and [14] show distinct spectral slope behaviors over 
the nucleus surface. We therefore propose three groups of ter¬ 
rains based on the strength of the spectral slope: 

- Low spectral slope, group 1. For this group, the spectral 
slope comprises between 11.0 and 14 %/(100 nm) at a phase 
angle of ~ 50° (Fig. HU) . it mostly corresponds to the blue 
color in Figs. [13] and [14] at different phase angles (10° and 
52°, respectively). For these regions, the albedo is generally 
higher. The regions Hathor, Hapi, and Seth are strongly cor¬ 
related with this group, although not the entire Seth region is 
blue. We might also add portions of the Ma’at (see Fig.[T3]at 
UT 17:26) and the Babi region (see Fig.[T4]at 23:19 UT). 

- Average spectral slope, group 2. For this group, the spectral 
slope comprises between 14 and 18%/(100 nm) at a phase 
angle of ~ 50° (Fig. [14]), it mostly corresponds to the cyan- 
green color in Figs. [13] and [14] Several regions are entirely 
correlated with this group (e.g., Anuket and Serqet), and oth¬ 
ers are partially correlated (e.g., Ma’at and Ash). 

- High spectral slope, group 3. For this group, the spectral 
slope is above 18%/(100 nm) at a phase angle of ~ 50° 
(Fig.m, it mostly corresponds to the orange and red colors 
in Figs. [ 13 ] and [14] Several regions of the nucleus are asso¬ 
ciated with this group, although only the Apis region is en¬ 
tirely correlated with group 3. In the small lobe, the regions 
around the Hatmehit depression have this high spectral slope 
(i.e.. Nut, Maftet, Ma’at, and Bastet), and a little of the inside 
of Hatmehit (Figure [14] at 01:19 UT). In the large lobe, the 
regions around the Imhotep depression are also correlated to 
the strong spectral slope (i.e.. Apis, Ash, and Khepry) and a 
little of the inside of Imhotep (Fig.[T4]at 06:19 UT). 

The three groups are observed in the two portions of the 
nucleus, the small and large lobes. As already stated (Sierks et 
al., 2015), there is no obvious color variability between the two 
lobes to support the contact-binary hypotheses, although at the 
same time, it is not ruled out either since two compositionally 
similar objects may have collided. Thus the color variations of 
the spectral slopes on the surface may be due to morphological 
and/or mineralogical properties. We explored this idea by rely¬ 
ing on the detailed analysis of the morphological properties of 
each region by El-Maarry et al. (2015). The group 1 regions are 
classified in various nomenclatures (i.e., brittle, dust cover, con¬ 
solidated), thus not highlighting a specific morphology. This is 
specifically true for the Hapi and Hathor regions, whose mor¬ 
phology is quite different. However, it is important to note that 
the Seth, Ma’at, Babi, Hathor, and Hapi regions are adjacent to 
each other, which concentrates the location of group 1 around 
the interface between the two lobes. These regions also have in 
common that activity has been spotted in most of them (Sierks 
et al., 2015, Vincent et al., 2015). Therefore, as shown for other 
comet nuclei (Sunshine et al., 2006; Li et al., 2013), the spectral 
slope of group 1 is most likely linked to water ice mixed with 
refractories on the surface, rather than a specific morphology of 
the nucleus. 

As previously noted, the albedo is generally correlated to the sur¬ 
face texture, but is not uniquely correlated to the spectral slope. 
In fact, smooth active regions such as Hapi and parts of Ma’at 
have both a very weak slope and high albedo (Figs. [8] and [14]). 
On the other side, the Imhotep region, which has a large portion 
of smooth material with higher albedo, did not show signs of ac¬ 
tivity during these observations, and its spectral slope belongs to 
groups 2 or 3. 
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Fig. 13. Spectral slope images for the observations on 1 August, at a phase angle of 9-10°. The slope is computed in the 535-882 nm range, after 
normalization at 535 nm, and it is in %/(100 nm). 


Spectral slope 



Fig. 14. Spectral slope images for the observations on 6 August, at a phase angle of ~ 50°. The slope is computed in the 535-882 nm range, after 
normalization at 535 nm, and it is %/(100 nm). 
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We also immediately noted that group 3 regions are generally 
located opposite the neck region. Both ends of the nucleus are 
characterized by very large depressions, whose origin might be 
different (see Auger et al. (2015) and Pajola et al. (2015) within 
this issue for more details). There is no evidence that these de¬ 
pressions are related to impacts (Pajola et al., 2015), thus the 
correlation of group 3 with ejecta deposits is not supported, also 
given the asymmetry of group 3 around these two depressions. 
Nonetheless, although only the Ma’at and Ash regions are clas¬ 
sified as "dust covered" (El-Maarry et al., 2015, Thomas et al., 
2015), Maftet, Apis, and Bastet are also adjacent to the depres¬ 
sion, and their surfaces appear rather smooth in the portions cor¬ 
related with group 3. Consequently, although the dusty regions 
are not clearly correlated to group 3, they are connected. This 
again emphasizes the relationship between surface texture and 
albedo. 

The original Philae landing site, Agilkia, is analyzed in detail 
in La Forgia et al. (2015), who reported a detailed geomorpho- 
logical and spectrophotometric study of this region and of the 
surrounding Hatmehit area from OSIRIS images acquired at a 
resolution of -0.5 m/px. These regions fall in the spectral groups 
2 and 3 according to our analysis. 

Another important aspect is that the groups’ spectral slope prop¬ 
erties are not correlated to a particular morphology that may ex¬ 
pose material from deeper inside the nucleus. In fact, the Aten 
depression belongs to group 2, which represents the average 
properties of the surface. The Nut depression, although it shows 
group 3 properties, is also significantly characterized by group 2 
properties. Therefore, the spectral slope variations of the surface 
do not show evidence of vertical diversity in the nucleus com¬ 
position, at least for the first tens of meters. However, it must be 
noted that sublimation lags, airfall, and surface transport com¬ 
plicate our interpretation, as these processes tend to homogenize 
the surface layer of the nucleus; below this uppermost layer, the 
subsurface might be different. 

5. Conclusions 

The data from comet 67P/Churyumov-Gerasimenko obtained 
from July to mid August 2014 with the OSIRIS imaging 
system show the nucleus with unprecedented spatial resolution. 
We have presented the results from the disk-averaged and 
disk-resolved analysis of the photometric properties and spec¬ 
trophotometry of the nucleus. Our main results are the following: 

- The phase function of the nucleus of 67P from disk-averaged 
reflectance in the 1.3°-54° phase angle range shows a strong 
opposition effect. The fit of the phase function with the HG 
model gives a G parameter of -0.13, implying a steep bright¬ 
ness dependence on the phase angle. The absolute magnitude 
reduced to the Bessel V filter is i7y(l, 1,0) = 15.74+0.02 
mag, and the linear slope (for a > 7°) is = 0.047+0.002 
mag/°, a value very similar to that found for comets Hartley 
2 and Tempel 1, and close to the average value for Jupiter- 
family comets. 

- We presented disk-integrated spectrophotometry in 20 NAG 
and WAG filters covering the 250-1000 nm range. The spec¬ 
tral behavior of the nucleus is red, meaning that the re¬ 
flectance increases with wavelength, and it is featureless in 
the visible and near-infrared range. This behavior is fully 
consistent with ground-based observations of the comet, and 
it is similar to those of other bare cometary nuclei. In the 
mid-UV region we see an increase of the flux that indicates 


a potential absorption band centered on - 290 nm, possibly 
due to SO 2 ice. However, this feature needs to be confirmed 
with other observations as the two WAG UV filters showing 
the flux enhancement are affected by pinhole defects. 

- The analysis of the spectral slopes versus phase angle, both 
from disk-averaged and disk-resolved data, shows a signifi¬ 
cant phase reddening for the nucleus, with the disk-averaged 
spectral slope increasing from 11%/(100 nm) to 16%/(100 
nm) in the 1.3°-54° phase angle range. This behavior may 
be attributed to a higher contribution from multiple-particle 
scattering at large phase angles and/or to surface roughness 
effects. The strong phase reddening suggests that the ef¬ 
fect of multiple scattering on the 67P nucleus may be non- 
negligible, even though the albedo is very low. 

- We modeled the disk-averaged reflectance in eight filters 
covering the 325-1000 nm range using the disk-integrated 
Hapke model formalism. We found no clear wavelength de¬ 
pendence of the g^, Bo and hg parameters, implying that the 
shadow-hiding effect must be the main cause of the oppo¬ 
sition surge. The geometric albedo at different wavelengths 
derived from Hapke modeling perfectly matches the comet 
spectrophotometry behavior. 

- Hapke (2002 and 2012) modeling on disk-resolved re¬ 
flectance at 649 nm gives parameters very close to those 
found for comets Wild 2 and Tempel 1 and is compatible 
with those found for comet Hartley 2, revealing that the 
photometric properties of these cometary nuclei are simi¬ 
lar. The geometric albedo derived from Hapke modeling is 
0.065+0.02 at 649 nm. Thus the 67P nucleus has a surface 
dark in absolute terms, but one of the brightest of the other 
cometary nuclei investigated by space missions (Li et al., 
2013). 

- The nucleus shows color and albedo variations across the 
surface: Hapi is ~ 16% brighter than the mean albedo over 
the surface, while the Apis and Seth regions are about 8-10% 
darker. 

- Based on local spectrophotometry, we see an enhancement 
of the flux at 701 and 743 nm that is clearly due to cometary 
emissions from the coma located between Rosetta and the 
nucleus. The potential sources of these emissions are H20^ 
in its (0,6,0) band at 700 nm, the (0,2,0)-(0,0,0) band at 743- 
745 nm, and the band system (0,6,0) of NH 2 , which overlaps 
with the water bands. 

- On the basis of the spectral slope and of the spectrophotome¬ 
try, we identified three different groups of regions, character¬ 
ized by a low, medium, and high spectral slope, respectively. 
The three groups are observed in the two lobes of the comets, 
and we did not note a significant color variability between the 
two lobes that would support the contact-binary hypothesis, 
although this does not exclude that the nucleus of 67P might 
have been formed by a collision of two objects in the past. 

- Generally, the spectral slope values are often anticorrelated 
with the reflectance, the brightest regions being also the bluer 
ones in terms of spectral slope. Hapi, in particular, is the 
brightest area in absolute terms, the main source of cometary 
activity at large heliocentric distances, and the spectrally 
bluer region. We interpret this behavior as due to a larger 
abundance of water ice in this region, which is often cast in 
shadow, not illuminated at perihelion due to seasonal effects, 
and thus a favored region to retain ices on its surface. 
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